Mapping the hydrodynamic response to the initial geometry in heavy-ion colhsions 
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We investigate how the initial geometry of a heavy-ion collision is transformed into final flow ob- 
servables by solving event- by- event ideal hydrodynamics with realistic fluctuating initial conditions. 
We study quantitatively to what extent anisotropic flow is determined by the initial eccentricity 
e„ for a set of realistic simulations, and we discuss which definition of e„ gives the best estimator of 
Vn- We find that the common practice of using an weight in the definition of e„ in general results 
in a poorer predictor of Vn than when using r" weight, for n > 2. We similarly study the importance 
of additional properties of the initial state. For example, we show that in order to correctly predict 
114 and V5 for non-central collisions, one must take into account nonlinear terms proportional to el 
and £2£3, respectively. We find that it makes no difference whether one calculates the eccentricities 
over a range of rapidity, or in a single slice at 2 = 0, nor is it important whether one uses an energy 
or entropy density weight. This knowledge will be important for making a more direct link between 
experimental observables and hydrodynamic initial conditions, the latter being poorly constrained 
at present. 

PACS numbers: 25.75.Ld, 24.10.Nz 
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I. INTRODUCTION 

Anisotropic flow [1] is one of the most important probes 
of ultrarelativistic nucleus-nucleus collisions. While early 
studies [2] focused on elliptic flow generated by the al- 
mond shape of the interaction region in non-central col- 
lisions, most of the recent activity concerns the effect of 
fluctuations in the initial geometry [3] . Such fluctuations 
result in fluctuations of elliptic flow [4], and also in new 
types of flow, such as triangular flow [5] and higher har- 
monics. These new flow observables have been recently 
measured at RHIC [6, 7] and LHC [8-11]. 

Flow phenomena are best modeled with ideal [12] or 
viscous [13] hydrodynamics. Event- by- event hydrody- 
namics [14] provides a natural way of studying flow fluc- 
tuations: one typically supplies a set of initial conditions, 
then evolves these initial conditions through ideal [14-18] 
or viscous [19] hydrodynamics, then computes particle 
emission at the end. Observables are finally averaged 
over a large number of initial conditions, much in the 
same way as they are averaged over events in an actual 
experiment. 

The largest source of uncertainty in these hydrody- 
namic models is the initial conditions [20, 21] — that is, 
the state of the system after which it has sufficiently ther- 
malized or isotropized for the hydrodynamic description 
to be valid. Several models of initial geometry fluctua- 
tions have been proposed [22-26] . The usual procedure is 
to choose one or two of these simple models for the initial 
conditions and calculate the resulting flow observables. 
Signiflcant progress has been made recently by simulta- 
neously comparing to several of the newly-measured flow 
observables. With this approach, hydrodynamic calcu- 
lations can be used to rule out a particular model of 
initial conditions if results do not match experimental 
data [6, 27, 28]. But it does not tell us why a particular 



model fails. In order to constrain the initial state directly 
from data, we need to identify which properties of the ini- 
tial state determine a given observable. These constraints 
can then provide valuable guidance in the construction 
of better, more sophisticated models of the early-time 
dynamics. 



It is well known that elliptic flow is largely determined 
by the participant eccentricity [4]. Teaney and Yan [29] 
have introduced a cumulant expansion of the initial den- 
sity profile, in which the participant eccentricity is only 
the first term in an infinite series, and they have sug- 
gested that the hydrodynamic response may be improved 
by adding higher-order terms, but to our knowledge their 
suggestion has never been checked quantitatively. Other 
expansions have also been suggested [30] . As for triangu- 
lar flow, ^3, symmetry considerations have been used to 
argue that it should be created by an initial triangularity 
£3, but several definitions of £3 are in use [5, 31] and it 
has never been investigated which is a better predictor 
of V3. Finally, it has been shown that higher harmon- 
ics [18, 32] V4 and W5 are in general not proportional to 
the corresponding £4 and £5. A possible better estimator 
was recently suggested [33], but it has not been checked 
quantitatively. 



The goal of this paper is to improve our understand- 
ing of the hydrodynamic response to initial fluctuations. 
We carry out event-by-event ideal hydrodynamic calcu- 
lations with realistic initial conditions and then quan- 
titatively compare the final values of w„ with estimates 
derived from the initial density profile. We are thus able 
to systematically determine the best estimators of fiow 
observables w„, n = 2-5 from the initial transverse den- 
sity profile. 
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II. A SYSTEMATIC APPROACH TO 
CHARACTERIZING THE HYDRODYNAMIC 
RESPONSE 

In hydrodynamics, the momentum distribution of par- 
ticles at tlic! end of the evolution is completely determined 
by initial conditions. Current models of initial conditions 
predict the system at early times to consist of flux tubes 
or other string-like structures that are extended longi- 
tudinally [34], with an approximate boost invariance at 
mid-rapidity. Most models also predict that the initial 
transverse flow, if any, is small [35, 36] (except for pos- 
sible fluctuations of the initial flow velocity [37]). Under 
these approximations, any observable is a completely de- 
terministic functional of the transverse energy density 
profile p{x,y). 

This density in turn can be completely characterized 
by a set of complex moments [29] 



p+q,p-<i 



j {x + iyY{x - iyyp{x, y)dxdy 
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Wo,o{r^'+''e'(f-«)'^}, (1) 



where {• • • } denotes an average value over the transverse 
plane weighted by p{x,y). Small values of the first in- 
dex (p + q) correspond to small powers of [fcj in a two- 
dimensional Fourier transform of p{x,y), and thus de- 
scribe large-scale structure, while moments with larger 
values are more sensitive to small-scale structure. The 
second index indicates the rotational symmetry of each 
moment. 

If the system has </> — >■ (j)+'jT symmetry, all odd moments 
(i.e., with p — q odd) vanish. In a symmetric heavy-ion 
collision, the symmetry between target and projectile im- 
plies an approximate <p ^ (p + n symmetry in a centered 
coordinate system defined by Wi^i = (used throughout 
this article). This symmetry is only broken by quantum 
fluctuations in the wavefunction of incoming nuclei [5]. 
Therefore, odd moments are typically small relative to 
even moments. Similarly, central heavy-ion collisions are 
rotationally symmetric except for fluctuations, so that 
all moments with p ^ q are small. For semi-central or 
peripheral collisions, however, the interaction area is al- 
mond shaped [2] , resulting in sizable moments in the 2nd 
Fourier harmonic p — q = 2. For example, the familiar 
participant eccentricity £2 and participant plane $2 are 
defined by [4] 



£26^'*= = - 



W2,2 _ {r^e^**} 



W2fi 



(2) 



The value of £2 is typically 0.3 in a semicentral heavy-ion 
collision: anisotropics arc small, and so it is natural to 
expect that the hydrodynamic response can be ordered 
into a Taylor's series. Higher-order even harmonics are 

smaller: The fourth-order anisotropy £4 = \W4^4\/W4fi is 
typically of order (£2)^ [38], so that it can be treated in 



practice as a higher order term in a Taylor series expan- 
sion. 

Any observable can generally be written as a function 
of moments of p. In this paper, we focus on anisotropic 
flow, Vn, which, along with the event-plane angle \E'„, is 
defined by 



(3) 



Here { • • • } denotes an average over the distribution of 
particle momenta in one event. In hydrodynamics, this 
is a smooth (boosted thermal) probability distribution. 
In a real- world collision, w„ and must be inferred from 
a finite sample of particles. The resulting statistical error 
makes an event-by-event determination of Vn impossible 
in current experiments — only event-averaged quantities 
are reliable. In theoretical calculations, however, we can 
accurately determine Vn and in every event in order 
to study precisely how they depend on initial conditions. 

The symmetries of u„ restrict what combinations of 
moments of the initial distribution it can depend on. 
For example, to first order in anisotropies, the rotational 
symmetry of u„ implies that it is a linear combination of 
moments in the same harmonic: 



— kn+2p,nWn+2p,n, 
p=0 



(4) 



where the coefficients kn+2p.n arc (dimensionful) func- 
tions of (the infinite set of) rotationally symmetric mo- 
ments W2m,o- The conventional eccentricity scaling of 
elliptic flow, n = 2, amounts to truncating the series to 
the first term, p = 0. I.e., the statement that V2 oc £2 
is a statement that the hydrodynamic response is sensi- 
tive only to the large-scale structure of the initial density 
distribution, with the response to small-scale structure 
damped in comparison. This statement has not been 
quantitatively tested until now. 

Teaney and Yan [29] have suggested that including a 
higher-order term p = 1 in addition to the lowest or- 
der may improve the accuracy (actually, they listed cu- 
mulants instead of moments, but it is equivalent). This 
conjecture will be checked quantitatively in Sec. IV. 

In this work we will use the following notation for the 
dimensionless eccentricity £„,« and the corresponding 
orientation angle ^m,n in a given event: 



(5) 



and we use the shorthand notations s„ 



^n,n- If m — n is even and positive, the numerator of 
Eq. (5) is Wm,n/Wofi- If is even, the denominator is 

Gubser and Yarom have proposed a different basis for 
the expansion [30], which can be seen as a partial re- 
summation of the infinite series (4). To first order in 

anisotropies, they write line™*" oc /nC™*"^, where /„ 
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and are solely determined by the initial density pro- 
file. The first harmonics are given by 



qre 



he-" 



1 + q-^r^ 
1 + q^r"^ 



(6) 



with q'^ = = T4^o,o/W^2,o- Strictly speaking, this 

expansion scheme is of obvious relevance only for a con- 
formal equation of state and for a particular initial den- 
sity profile falling more slowly at large r than realistic 
profiles. Nevertheless, we find it instructive to test how 
this expansion compares with conventional eccentricity 
scaling with realistic initial conditions. 

Finally, one can also consider moments of the entropy 
density profile instead of energy density, which we also 
test in the following. 

The goal here is to determine which moment or com- 
bination of moments serve as a best estimator for flow 
observables Vn with n = 2-5. We shall see in particular 
that V4 and are not well described by the leading-order 
expansion (4) and require nonlinear terms, which are also 
constrained by symmetry. Note that symmetry consider- 
ations alone exclude a linear mixing between the second 
and third harmonics, as proposed in Ref. [39] . 



III. DETERMINING THE BEST ESTIMATOR 

OF Vn 

The goal of this work is to test to what extent u„ and 
^'n are correlated with quantities derived from the initial 
transverse density distribution, such as e„ and <&„. Pre- 
viously, the correlation of anisotropic flow with the initial 
geometry has been studied by plotting the distribution 
of — $„ [16, 18, 31] or by displaying a scatter plot of 
Vn versus e„ [40]. In this paper, we carry out a global 
analysis which studies both aspects simultaneously and 
quantitatively. 

For a given event, we write 



Vr,.e 



(7) 



where k is an unknown proportionality constant. The 
first term in the right-hand side defines the estimate for 
Vn from the initial eccentricity, and the last term £ is the 
difference between the calculated flow and the proposed 
estimator, or the error in the estimate (note that £ is 
complex). No known estimator can perfectly predict Vn 
in every event (for example, two events with the same 
triangularity can be constructed to have very different 
triangular flow [32, 41]). The best estimator, then, should 
be deflned as the one that minimizes the mean square 
error (|f P), where (• • • ) denotes an average over events 



in a centrality class. A straightforward calculation shows 
that the best value of k is 



k = 



{SnVn COs(n(^„ - $„))) 
(4) 



(8) 



Inserting Eq. (8) into Eq. (7), one obtains the best esti- 
mator of Vn from En in a centrality class. Using Eqs. (7) 
and (8), one finally derives the mean-square error 



{\£\') = (vl) - e (el) . 



(9) 



This sho ws that the rms value of the best estimator, 
\k\y^ (e^j), is always smaller than the rms value of w„. 
In the next section, we compute k^J{e^ / yjjv^ for var- 
ious definitions of e„ and several values of n. The closer 
the ratio to 1, the better the estimate. Using Eq. (9), 
a ratio of 0.95 corresponds to a rms error of 31%. A 
change of sign in the ratio signals that the estimator is 
anticorrelated to Vn- 

According to the discussion in Sec. II, an improved esti- 
mator may be obtained by adding more terms in Eq. (7), 
e-g-, 



Vr,e 



(10) 



where e'n and <I>'j are other quantities determined from 
the initial density profile (for example, the next higher 
cumulant). The best estimator is now given by the fol- 
lowing system of equations 

{snVn cos(n(*„ - $„))) = A;'(e„4 cos(n($J, - $„))) 

+k{el) 

(e>„cos(n(*„ - $;))) = fc(£„£' cos(n($„ - ^'J)) 

+k'{e'b, (11) 

which can be solved for k and k' . The rms value of the 
best estimator and the rms error are related by an equa- 
tion analogous to Eq. (9): 



(12) 



One can show that the rms error is always smaller with 
two terms, Eq. (12), than with only one of the terms, 
Eq. (9). This is intuitive if one thinks of Eqs. (7) and 
(10) as fits to Vn- adding more parameters improves the 
quality of the fit. 



IV. RESULTS 

We simulate Au-Au collisions at the top RHIC energy 
using the hydrodynamic code NeXSPheRIO [14]. NeX- 
SPheRIO solves the equations of relativistic ideal hydro- 
dynamics using initial conditions provided by the event 
generator NeXus [42]. Fluctuations in initial conditions 
are studied by generating 150 NeXus events in each of the 
10% centrality classes studied, and solving the equations 
of ideal hydrodynamics independently for each event. In 
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addition, 115 NeXus events with zero impact parame- 
ter were used in order to study very central collisions. 
NeXSPheRIO provides a good description of rapidity 
and transverse momentum spectra [43], and elliptic flow 
V2 [44]. In addition, it reproduces the long-range struc- 
tures observed in two-particle correlations [45] . 

The code NeXSPheRIO emits particles at the end of 
the hydrodynamical evolution using a Monte-Carlo gen- 
erator. Anisotropic flow w„, and the corresponding event- 
plane angle are defined from Eq. (3), where {• • • } de- 
notes an average over all particles in the pseudorapidity 
interval — 1 < r/ < 1. 

This work requires an accurate determination of w„ in 
each hydrodynamic event, and so associated with each 
initial condition, we generate approximately 6 x 10^ par- 
ticles by computing particle production with N Monte- 
Carlos. This allows for a much better event plane resolu- 
tion and much smaller statistical error (the actual mul- 
tiplicity in an event is « 6 x 10^/iV). We compute V2 
to V5. We do not compute vi because it changes sign 
as a function of transverse momentum [40], so that an 
average with equal weighting such as in Eq. (3) is not 
appropriate [46]. An analysis of directed flow is left to 
future work. The relative statistical errors on V2 to U5 
in a given event are 3.7%, 5.7%, 9.8%, 20%, respectively. 
The rms error on the event planes ^2 to ^5 are 1°, 1°, 
1.5°, 2°. This means that the event-plane resolution [47] 
for is as large as 0.98, and even closer to 1 for all other 
harmonics. 

In such a 3-dimensional calculation, there is more than 
one way to define the transverse energy density profile 
that is used as a weight when calculating the eccentric- 
ities in Eq. (5). We show results obtained by averaging 
over the transverse energy density profile at z = (i.e., 
central space-time rapidity 77s = 0). Though the results 
are not shown, we have found that averaging over the 
space-time rapidity interval — 1 < t^s < 1 results in pre- 
dictors of equal quality. 

It should be noted that the error from these predictors 
is likely to be larger in our calculations than in many 
others, for several reasons. First, our hydrodynamical 
calculations are based on NeXus initial conditions which 
contain (fluctuating) initial flow, as well as longitudinal 
fluctuations, and so the final flow measured in a partic- 
ular pseudorapidity window is not entirely determined 
by the initial transverse geometry. In addition, there are 
statistical fluctuations from the finite number of particles 
generated at the end of the hydro evolution. These issues 
set a limit on the rms error introduced in Sec. Ill, which 
cannot go to zero. A hydrodynamic calculation with less 
or no initial flow, in 2+1 dimensions, or that calculates 
flow from a continuous distribution at freeze out, though 
perhaps less realistic, will likely result in a smaller error 
for the same estimator. Likewise, a non-zero viscosity 
may cause higher-order cumulants to decrease in impor- 
tance, improving the predictive power of the lowest mo- 
ments En- In this sense, these results represent something 
of a worst-case scenario. 
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FIG. 1. (Color online) Best estimator for elliptic flow divided 
by rms V2 for various combinations of moments of the ini- 
tial energy density. The leftmost points correspond to 115 
events with exactly zero impact parameter. Error bars repre- 
sent statistical uncertainty from the finite number of events. 
Diamonds (squares) have been given an x-offset of 1.5 (-1.5) 
for readability. 



A. Elliptic flow 

Elliptic flow is usually thought of as a hydrodynamic 
response to the initial eccentricity: 1)26^'*^ = A:e2e^**^, 
where the resulting event plane "^2 approximately coin- 
cides with the participant plane $2 [4]. The right-hand 
side of this equation defines an estimator of V2- 

Fig. 1 displays the ratio of the rms value of the best 
estimator to the rms value of the calculated V2 as a func- 
tion of centrality. The proportionality constant k is de- 
termined using Eq. (8) independently for each centrality 
class. For central collisions, where all anisotropics are due 
to fluctuations, the best estimator is able to reproduce 
over 80% of V2 ■ This means that the participant eccen- 
tricity correctly captures the physics of V2 fluctuations, 
but that the event-plane ^'2 fluctuates around the par- 
ticipant plane $2 [16, 18, 31], and/or that V2 has sizable 
fluctuations for a given £2- For mid-central collisions, el- 
liptic flow is driven by the almond-shaped overlap area, 
therefore fluctuations are smaller and the estimate is bet- 
ter. The value of fc slightly decreases with centrality, but 
very slowly. It is approximately equal to 0.16. Note that 
k does not represent the ratio of the magnitudes V2/£2 
in a typical collision, and should not be compared with 
V2/ £2 obtained from smooth initial conditions, nor to the 
ratio of the average elliptic flow to the average eccentric- 
ity V2/S2 which are both larger [27]. 

As explained in Sec. II, the participant eccentricity 
is one term out of an infinite series of moments (or 
cumulants) allowed by symmetry. There is no funda- 
mental reason why the first term in Eq. (4) must be 
more important than higher-order terms. In order to 
check quantitatively this issue, we define another estima- 
tor of V2, corresponding to the term p = 1 in Eq. (4): 
f2e^'*^ = fc£4.2e^**'''^ . The difference with usual partic- 
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ipant eccentricity scaling is that larger values of r are 
given more weight. As shown in Fig. 1, this estimate is 
essentially as good as the usual participant eccentricity. 
A closer look reveals that it is slightly better for central 
collisions, and slightly worse for non-central collisions. 
This means that V2 is driven more by the periphery of 
the fireball for central collisions than for peripheral colli- 
sions. The result that, for central collisions, anisotropic 
flow Vn is sensitive to the geometry of the outer layers 
of the system is in agreement with Ref. [32, 48]. There 
is a limit, however. We have checked that £5,2 gives a 
worse estimate than £4^2 for all centralities, indicating 
that higher-order moments, and thus the extreme pe- 
riphery, are indeed less important. 

Using entropy density instead of energy density to cal- 
culate moments also gives a somewhat greater weight to 
larger r. Indeed we have found that, calculating £2 with 
an entropy density weight gives results (not shown) that 
are between the results for £2 and £4^2 — that is, slightly 
better for central collisions and slightly worse for periph- 
eral collisions, but in general the result is very close to 
the result for £2 calculated using energy density. In gen- 
eral either predictor appears to be essentially as good as 
the other. 

Next, we test if the quality of the estimator is improved 
by combining £2 and £4,27 as in Eq. (10). The improve- 
ment is marginal, which means that adding the next term 
in the cumulant expansion [29] does not significantly im- 
prove the determination of the event plane from initial 
conditions. This implies that the small-scale structure 
of the initial conditions is unlikely to be responsible for 
the part of elliptic flow that is not explained by £2, and 
that most of the deviation from being a perfect predictor 
is likely coming from another source — e.g., nonlinear 
terms or fluctuating initial flow. 

Finally, we test Gubser's estimator, 2nd line of Eq. (6). 
This particular quantity gives less weight to the periph- 
ery. This makes the estimator much worse for central 
collisions, as expected from the discussion above, but it 
is also worse at all other centralities. 

Overall, £2 calculated with energy or entropy density 
weighting is a very good predictor of V2, while other quan- 
tities that are significantly more sensitive to the periph- 
ery compared to the center of the system, or vice versa, 
are worse. 



B. Triangular flow 

Similar to elliptic fiow and eccentricity, triangular flow 
V3 is thought to be a hydrodynamic response to an initial 
"triangularity" in the initial state: to deflne the trian- 
gularity, Alver and Roland [5] originally suggested £2,3 
(following the notation of Eq. (5)). However, the nu- 
merator cannot be simply expressed in terms of the mo- 
ments (1). Recently it has been more common to use 
£3,3 [31] J which we denote by £3. In this case, the nu- 
merator is Ws.s/Wq^o, but the denominator {r^} is not 
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FIG. 2. (Color online) Best estimator for triangular flow di- 
vided by rms 113 for various definitions of £3. The leftmost 
points correspond to 115 events with exactly zero impact pa- 
rameter. Error bars represent statistical uncertainty from the 
finite number of events. Triangles (squares) have been given 
an x-offset of 1.5 (-1.5) for readability. 



a simple moment. One could instead replace {r^} with 
a power of the lowest moment {r^}^^^. This gives an al- 
most indistinguishable result in our analysis, so we only 
show the curve for the "standard" denominator.^ 

Figure 2 shows the ratio of the rms value of the best es- 
timator to the rms value of the calculated V3. The trian- 
gularity £2.3 is a worse predictor of than £3 below 30% 
centrality, but slightly better above 40%. As with V2, the 
second-lowest moment £5,3 is a slightly better predictor 
for central collisions, but worse for non-central collisions, 
while £7^3 (not shown) is worse everywhere. This again 
signals a somewhat stronger sensitivity to the periphery 
of the collision region in central collisions than in periph- 
eral collisions, though again the moment fs, which has a 
strong sensitivity to the interior, is worse at all centrali- 
ties. 

As with elliptic flow, replacing an energy density 
weight with an entropy density weight in the calculation 
of £3 gives results (not shown) that are slightly better for 
central collisions and worse for peripheral collisions, but 
that are essentially equivalent. 

Finally, using a sum of the lowest two moments £3 
and £5^3 (not shown) reproduces the highest points on 
the figure. I.e., it shows no improvement over the term 
that is individually the best predictor, except above 40% 
centrality, where it follows the £2,3 result. 

Thus, £3 is a very good predictor, with slightly too 
little sensitivity to the periphery in central collisions and 
slightly too much in peripheral collisions, but quantities 
with too-different r-dependence are everywhere worse. 



^ Values of £3 are slightly larger with {r^y^/^ than with {r^} de- 
nominator, and £3 is no longer bounded by 1, but this is almost 
exactly compensated by the smaller value of k from Eq. (8). 



6 



-•-£4 

^6.4 
-A- ^2,4 




% centrality 



FIG. 3. (Color online) Best estimator for quadrangular flow 
divided by rms 114 for various choices of the estimator. Nega- 
tive values indicate an anticorrelation with the estimator. The 
leftmost points correspond to 115 events with exactly zero im- 
pact parameter. Error bars represent statistical uncertainty 
from the finite number of events. Triangles (squares) have 
been given an x-offset of 1.5 (-1.5) for readability. 



C. Quadrangular flow 

Like the triangularity, the quadrangularity has incon- 
sistent definition in the literature. Some authors use 
£2,4 [27], others use £4^4 (which we denote by £4) [38]. 
Another possible choice is, in the latter case, to replace 
{r^} with {r^}^ in the denominator. In this case, the es- 
timator is equally good, and so the results are not shown 
with the rest of the results in Fig. 3. 

These results differ qualitatively from our results for 
V2 and U3. For V2 and ^3, all the estimators we have 
tested give good results, for all centralities. By contrast, 
£4 gives reasonable predictions only for central collisions. 
The agreement becomes much worse for peripheral colli- 
sions, in agreement with previous analysis [18]. In fact, U4 
is anticorrelated with £4 for the most peripheral points. 
Despite suggested scaling similarities between the cen- 
trality dependence of £4 and U4 [38], this shows that £4 
cannot be used as an estimator of V4 on an event-by-event 
basis for non-central collisions. 

When using an weight (£2,4), the result is signifi- 
cantly worse for central collisions. Using the next highest 
moment £5,4 is slightly better for central collisions, but 
all higher moments are worse. 

For peripheral collisions, the asymmetry of the nuclear 
overlap region, causes £2 to be significantly larger than 
other moments such as £4. This raises the possibility that 
non-linear terms involving £2 may be important. The 
first such term allowed by symmetry is proportional to 
el — that is, k{e2e'^"^^y [33]. Fig. 3 shows that, in- 
deed, this term alone provides a reasonable estimator for 
non-central collisions. More interestingly, including both 
terms, i.e., 

V4e^''^^ = ke^e^"^^ + k'ele^"^' , (13) 
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FIG. 4. (Color online) Best estimator for pentagonal flow 
divided by rms 115 for various choices of the estimator. Nega- 
tive values indicate an anticorrelation with the estimator. The 
leftmost points correspond to 115 events with exactly zero im- 
pact parameter. Error bars represent statistical uncertainty 
from the finite number of events. Triangles (squares) have 
been given an x-offset of 1.5 (-1.5) for readability. 

results in an excellent predictor for all centralities. For 
central collisions, £4 and £2 are both small and of the 
same order of magnitude, so that the first term domi- 
nates. For all other centralities, both terms are of com- 
parable magnitudes, and the combination gives a much 
better result than either individual term. 

Note that the rms W4 was first measured in 2011 [6, 9- 
11]. In earlier analyses, V4 was determined with respect 
to the event plane from elliptic flow [49, 50]. The mea- 
sured quantity is then (z;4e"***''w|e"^'*2)/(„2^ [51]^ ^ot 
the rms v^. 



D. Pentagonal flow 

Figure 4 shows the results for v^, which are very sim- 
ilar to the results for V4. Here, the non- linear term £2£3 
becomes important even for central collisions. The pre- 
dictor 

„5e5**^ = feseS**^ + /fc'£2e2'*= ese^"^' , (14) 
does an excellent job at all centralities. 



V. CONCLUSIONS 

In this work, we have quantitatively tested to what 
extent anisotropic flow can be predicted from the ini- 
tial density proflle in event-by-event ideal hydrodynamics 
with realistic initial conditions. We have shown that the 
participant eccentricity £2 gives a very good prediction of 
elliptic flow for all centralities. We have also shown that 
the definition of £3 with weights [29, 31] gives a better 
prediction of triangular fiow than the previous definition 
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with weights. Gubscr's moments [30] give worse re- 
sults for both V2 and . Higher harmonics W4 and W5 can 
be well predicted from the corresponding eccentricities 
£4 and £5 (again defined with r"^ and weights rather 
than with weights) only for central collisions. For 
noncentral collisions, a good predictor of V4 must include 
two terms, proportional to £4 and £3. Likewise, ^5 has 
contributions proportional to £5 and £2£3. Defining the 
eccentricities with energy or entropy density, or using the 
density at a midrapidity slice or over a finite longitudi- 
nal range is largely a matter of preference, and does not 
make a significant difference. 

These results provide an improved understanding of 
the hydrodynamic response to the initial state in realis- 
tic heavy-ion collisions, and provide a more direct link 
between experimental data and properties of the initial 



stage of the collision. This will allow for the construction 
of more realistic models for the early-time collision dy- 
namics, and thus a significant reduction in the systematic 
uncertainties of extracted bulk properties of the system. 
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